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1 Introduction 



The recent data of Super-Kamiokande (SK) have given a new impetus to the 
atmospheric neutrino problem and a possible interpretation in terms of neutrino 
oscillation. Moreover the high statistics of SK makes it possible to study the zenith- 
angle dependence of the neutrino flux from which one can conclude that the i/^'s 
show signs of oscillation but the v^. events are consistent with the no-oscillation 
hypothesis. Independently the results from the reactor experiment CHOOZ dis- 
favours the z/^ — oscillation hypothesis . On the other hand large angle z/^ — 
or — Vs being a sterile neutrino) solution continues to give a good fit to the 
data. Nevertheless effort has been on to try out other possibilities to explain the 
anomaly observed in SK and one among these is neutrino decay p, |]. In it was 
shown that neutrino decay gives a poor fit to the data. However they considered 
neutrinos with zero mixing. Barger et al. considered the situation of neutrino de- 
cay in the general case of neutrinos with non-zero mixing angle They showed 
that the neutrino decay fits the L/E distribution of the SK data well. The Am^ 
taken by them was > 0.1 eV"^ so that the Am^ dependent term averages out. As 
pointed out in such a constraint on Am^ is valid when the unstable state decays 
into some other state with which it mixes. If however the unstable state decays 
into a sterile state with which it does not mix then there is no reason to assume 
Am^ > 0.1 e1/2. 

In this paper we present our results of two-flavour z/^ — Ur oscillation and neu- 
trino decay solutions to the atmospheric neutrino problem by doing x^-fit to the 
848 days of sub-GeV and multi-GeV Super-Kamiokande data p. We have also 
presented the results of x^-fit to the 535 days SK data and have compared it with 
the results for the new data. For the neutrino decay analysis we take the most 
general case of neutrinos with non-zero mixing and consider two pictures 

• Am^ > 0.1 eV"^ (scenario (a)) 

• Am^ unconstrained (scenario (b)) 

We also explicitly demonstrate the behavior of the up-down asymmetry parameters 
P, 1^ in both scenarios. 

Our analysis shows that scenario (a) is ruled out at 100%(99.99%) C.L. by 
the 848(535) days of SK data. However if we remove the constraint on Am^ and 
consider the possibility of decay into a sterile state then one can get an acceptable 
flt for Am^ ~ 0.001e1/2 ^nd sin^ 29 large. 

The plan of the paper is as follows. In section 2 we present our results for 
two-generation — oscillation analysis. In section 3.1 we present our results for 
the neutrino decay solution constraining Am^ to be > O.leV^^. In section 3.2 we 
do a three parameter analysis by removing the constraint on Am?. In section 
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4 we perform a comparative study of the three cases and indicate how one can 
distinguish experimentally between the scenario (b) and the z/^ — v^- oscillation 
case though both give almost identical zenith-angle distribution. 



2 i^^ — i^T oscillation 

In the two-flavour picture the probability that an initial z/; of energy E remains a 
ui after traveling a distance L in vacuum is 

P,,,, = 1 - sin ^29 sm\7iL/Xosc) (1) 

where 9 is the mixing angle between the two neutrino states in vacuum and Xosc 
is the oscillation wavelength deflned as, 

A„,,= (2.5M (2) 
GeV Am^ 

where Am^ denotes the mass squared difference between the two mass eigenstates. 
The expected number of / (e or fi) like 1 ring events recorded in the detector in 
presence of oscillations is given by 

max r + l /• + ! 1 /•27r 

Ni = riT dE I dEi / dcosip / dcosC, — / d(f) 

Jo "'(-B,)mm J-i J-i 2tt Jo 

d'Fi{E,0 d'aiiE, El, cos ^P) 



dE dcosC, dEi dcosip 



-eiEi). P.,.,iE,0. (3) 



riT denotes the number of target nucleons, E is the neutrino energy, Ei is the 
energy of the final charged lepton, is the angle between the incoming neutrino 
ui and the scattered lepton ^ is the zenith angle of the neutrino and (p is the 
azimuthal angle corresponding to the incident neutrino direction (the azimuthal 
angle relative to the ip has been integrated out). The zenith angle of the charged 
lepton is then given by 

cos 9 = cos C, cos ip + sin^ cos (j) sin ip (4) 

d'^ Fi /dEd cos is the differential flux of atmospheric neutrinos of type ui, d'^ai/dEid cos ip 
is the differential cross section for uiN —>■ IX scattering and e{Ei) is the detection 
efficiency for the 1 ring events in the detector. The efficiencies that were available 
to us are not the detection efficiencies of the charged leptons but some function 
which we call e{E) defined as P] 

<E) = 'fl'L (5) 
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Puiui is the survival probability of a neutrino flavour / after traveling a distance L 
given by, 



L = ^{Re + hf - Rj?<iTi^i- Recosi (6) 

Re being the radius of the earth and h is the height of the atmosphere where the 
neutrinos are produced. We use the atmospheric neutrino fluxes from 0. For the 
sub-GeV events the dominant process is the charged current quasi-elastic scattering 
from free or bound nucleons. We use the cross- sect ions given in |]T0[. The events 



in multi-GeV range have contributions coming from quasi-elastic scattering, single 
pion production and multi pion production and we have used the cross-sections 
given in |ll|]. For the multi-GeV events we assume that the lepton direction G 



is the same as the incoming neutrino direction ^. But actually they are slightly 
different. We simulate this difference in the zenith angles by smearing the angular 
distribution of the number of events with a Gaussian distribution having a one 
sigma width of 15° for /i type events and 25° for the e type events [|12|. For 
the sub-GeV events, difference in direction between the charged lepton and the 
neutrinos are exactly taken care of according to eq. (|^) and (|[). 

To reduce the uncertainty in the absolute flux values the atmospheric neutrino 
measurements are usually presented in terms of the double ratio 

where MC denotes the Monte-Carlo simulated ratio. Different calculations agree 
to within better than 5% on the magnitude of this quantity. We use a similar 
quantity R, where 

WW)\ — ■ 

V-''^/-''e7| no— osc 

The quantities N^,^ are the numbers of e-like and /i-like events, as per eq.(|]). 
The numerator denotes numbers obtained from eq.@, while the denominator the 
numbers expected with the survival probability as 1. 

At the detector, the neutrino flux come from all directions. Thus, the total 
path length between the production point in the atmosphere and the detector 
varies from about 10 km to 13,000 km depending on the zenith angle. Neutrinos 
with zenith angle less than 90° {'downward neutrinos') travel a distance of ~ 10 
- 100 km from their production point in the atmosphere to the detector while the 
neutrinos with larger zenith angles {'upward neutrinos') cross a distance of up to 
~ 13,000 km to reach the detector. Apart from altering the flavour-content of 
the atmospheric neutrino flux, oscillations could lead to the following effect: if the 
oscillation length is much longer than the height of the atmosphere but smaller 
than the diameter of the earth, only upward neutrinos coming from the opposite 
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side of the earth will have significant oscillations. These would show up as an up- 
down asymmetry in the event distribution. SK has enough statistics to study these 
up-down flux asymmetry. They divide the (—1, +1) interval in cosG in five equal 
bins: (-1.0,-0.6), (-0.6,-0.2), (-0.2, +0.2), (+0.2, +0.6), (+0.6, +1.0) and give 
the number of events in each bin. The first two bins correspond to the upward 
neutrinos and the last two bins correspond to the downward neutrinos. To probe 
the up-down flux asymmetries we use the parameter Y defined in 0, 

Here Nf^''^ denotes the number of /-type events produced in the detector with 
zenith angle cosO < —0.2, i.e. the upward neutrino events while Nj^^''^ denotes 
the number of /-type events for cosB > 0.2 i.e. events coming from downward 
neutrinos. The central bin has contributions from both upward and downward 
neutrinos and is not useful for studying the up-down asymmetry. 
We minimize the function defined as [^] 




(10) 



where the sum is over the sub-GeV and multi-GeV cases. The experimentally ob- 
served rates are denoted by the superscript "exp" and the theoretical predictions 
for the quantities are labeled by "th" . AR^^^ is the error in R obtained by combin- 
ing the statistical and systematic errors in quadrature. AV^^p corresponds to the 
error in Y. For this we take only the statistical errors since these are much larger 
compared to the systematic errors. We include both the e-like and the /x-like up- 
down asymmetries in the fit so that we have 4 degrees of freedom (6 experimental 
data - 2 parameters) for the oscillation analysis in the two parameters Am^ and 
sin^ 29. 

The use of these type of ratios for the analysis test has been questioned 
in ||13[ because the error distribution of these ratios is non-Gaussian in nature. 
The alternative is to use the absolute number of e or type events taking into 
account the errors and their correlations properly H]. However as has been 



shown in the use of the i?'s and F's as defined above is justified within the 
3(7 region around the best-fit point for a high statistics experiment like SK and 
provides an alternative way of doing the x^-analysis. A comparison of the results 
of with those obtained in []T2|, [T^ shows that the best-fit points and the allowed 
regions obtained do not differ significantly in the two approaches of data fitting. 
The advantage of using the ratios is that they are relatively insensitive to the 
uncertainties in the neutrino fluxes and cross-sections as the overall normalization 
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factor gets canceled out in the ratio. We have included the Ye in our analysis 
because to justify the — Ur oscillation scenario, it is necessary to check that 
including the data on electron events gives a low value and hence it is the standard 
practice to include these in the x^-analysis 0, |12|, [1^. The data that we have used 



are shown in Table 1 which corresponds to the 848 days ^ and the 535 days [0 
of data. 

Table 1: The SK data used in this analysis. 



Quantity 


848 days data 


535 days data 


Sub - GeV 


Multi - GeV 


Sub - GeV 


Multi - GeV 




0.69 


0.68 


0.63 


0.65 




0.05 


0.09 


0.06 


0.09 


•yexp 


0.74 


0.53 


0.76 


0.55 


/\yexp 


0.04 


0.05 


0.05 


0.06 


e 


1.03 


0.95 


1.14 


0.91 


A yexp 
e 


0.06 


0.11 


0.08 


0.13 



For the 2 flavour z/^ — u^- oscillation the Xmin ^^at we get is 1.21 with the 
best-flt values as Am^ = 0.003 eV"^ and sin^ 26 = 1.0. This provides a good flt 
to the data being allowed at 87.64% C.L. If we use the 535 days data then the 
Xmin th^t S^^ is 4.25 with the best-flt values as Am^ = 0.005 eV"^ and sin^ 29 = 
1.0, the g.o.f being 37.32%. Thus the flt becomes much better with the 848 days 
data with no signiflcant change in the best-flt values. Though we have used a 
different procedure of data fltting, our results agree well with that obtained by the 
SK collaboration]]. In flg. 1 we show the 90% C.L. (x^ < Xmin + 4.61) and the 
99% C.L. (x^ < Xmin + 9-21) allowed region in the (Am^ , sin^ 26 ) plane for the 
~ oscillation hypothesis using the latest SK data. 



3 Neutrino decay 

The neutrino decay hypothesis assumes that there is an unstable component in 
Vfj, (say V2) which decays into one of the lighter states (say z/3). Experimental 
considerations constrain Vf. to decouple from V2 and it's decay partners, so that 

Ve ~ vi (11) 

i/^ ^ z/2 cos 6* + z/3 sin ^ (12) 

^The best-fit values that the SK coUaboration has got for the 848 days data are H] Am^ ~ 
0.003 eV^, sin^ 20 = 0.995 and x^/d-o.J = 55.4/67. This corresponds to a g.o.f of 84.33%. 
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From (|T2|) the survival probability of the z/^ of energy E, with an unstable compo- 
nent 1^2, after traveling a distance L is given by, 

P^^^^ = sin^e + cos^eexp(-47rL/Ad) 

+ 2 sin^ 9 cos 6'exp(— 27rL/A(i) cos(27rL/Aosc) , (13) 

where is the decay length (analogous to the oscillation wavelength given by eq. 
d^)) defined as, 

Ad = 2.5km-^^ (14) 
Gel/ a 

and a = ^2/70, m2 being the mass of the state z/2 and tq the decay lifetime. The 
Xosc appearing in eq. (0) is the wavelength of oscillations as defined in eq. (0) 
with Am^ = m| — m|. 

3.1 Am^ > O.leV^ 

If the unstable component in the z/^ state decays to some other state with which 



it mixes then bounds from K decays imply Am > O.leV [16]. In this case the 



cos(27rL/Aosc) term averages to zero and the probability becomes 

P,^,^ = sin^ 9 + cos^ ^exp(-47rL/Ad) . (15) 

In figs. 2 and 3 we show the variation of R and Y with a for various values 
of sin^ 9 for the sub-GeV and multi-GeV cases. For higher values of a, the decay 
length Ad given by eq. (0) is low and the exponential term in the survival proba- 
bility is less implying that more number of neutrinos decay and hence R is low. As 
a decreases the decay length increases and the number of decaying neutrinos de- 
creases, increasing R. For very low values of a the exponential term goes to 1, the 
neutrinos do not get the time to decay so that the probability becomes 1 — | sin^ 29 
and remains constant thereafter for all lower values of a. This is to be contrasted 
with the z/^ — u^- oscillation case where in the no oscillation limit the sin^(7rL/Aosc) 
term and the survival probability — > 1. For multi-GeV neutrinos since the 
energy is higher the A^ is higher and the no decay limit is reached for a larger value 
of a as compared to the sub-GeV case. This explains why the multi-GeV curves 
become flatter at a higher a. The behavior of the up-down asymmetry parame- 
ter is also completely different from the only oscillation case |]T^. In particular 
the plateau obtained for a range of Am^ which was considered as a characteristic 
prediction for up-down asymmetries is missing here. For the decay case even for 
a as high as 0.001 eV"^, the decay length A^ = 2500 (E/GeV) km so that the 
exponential term is 1, there is almost no decay for the downward neutrinos and 
the survival probability is P = 1 — | sin^ 29 while the upward going neutrinos have 
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some decay and so Y is less than 1. As a decreases, the increases, and the 
fraction of upward going neutrinos decaying decreases and this increases Y. For 
very small values of a even the upward neutrinos do not decay and Y ^ 1 being 
independent of 9. 

We also perform a analysis of the data calculating the "th" quantities in 
([T0|) for this scenario. The best-fit values that we get are a = 0.33 x 10~^ in eV"^ 
and sin^ 6* = 0.03 with a Xmm 49.16. For 4 degrees of freedom this solution is 
ruled out at 100% C.L. The best-fit values for the 535 days of data that we get 
are a = 0.28 x lO^'* in eV^ and sin^ 6 = 0.08 with a Xmm of 31.71. For 4 degrees 
of freedom this solution is ruled out at 99.99% C.L. |]18|. Thus the fit becomes 
worse with the 848 days data as compared to the 535 days data. We have marked 
the R and Y corresponding to the best-fit value of the parameters a and sin^ 9 
in figs. 2 and 3. It can be seen that the best-fit value of R for the sub-GeV 
neutrinos is just below and that for the multi-GeV neutrinos is just above the ±lcr 
allowed band of the SK 848 days of data. The up-down asymmetry parameter 
Y is quite low for the sub-GeV neutrinos and extremely high for the multi-GeV 
neutrinos as compared to that allowed by the data. The fig. 2 shows that for the 
sub-GeV neutrinos the data demands a lower value of a while from fig. 3 we see 
that the multi-GeV neutrinos need a much higher a to explain the SK data. It 
is not possible to get an a that can satisfy both the sub-GeV and the multi-GeV 
SK data, particularly it's zenith angle distribution. In this scenario, decay for the 
sub-GeV upward neutrinos is more than that for the multi-GeV upward neutrinos 
(downward neutrinos do not decay much) and as a result Y for sub-GeV is lower 
than the Y for multi-GeV, a fact not supported by the data. Since the 848 days 
data needs even lesser depletion of the sub-GeV flux as compared to the multi-GeV 
flux, the fit gets worse. 



3.2 Am^ unconstrained 

In this section we present the results of our x^-analysis removing the constraint 
on Am^. This case corresponds to the unstable neutrino state decaying to some 
sterile state with which it does not mix The probability will be still given by 
eq. dg). 

In fig. 4 and 5 we plot the R vs. Am^ and Y vs. Am^ for the sub-GeV and 
multi-GeV data for a = 0.3 xlO~^eV^^ (which is the best-fit value we get for the 
848 days data) and compare with the curve obtained for the best-fit value of sin^ 9 
(=0.5) for the only oscillation case (solid line). For the best-fit value of a that we 
get, the downward neutrinos do not have time to decay while the upward neutrinos 
undergo very little decay. Thus the curves are very similar in nature to the only 
oscillation curves. In the sub-GeV case (fig. 4), for high values of Am^ around 
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0.1 eV'^ both upward and downward neutrinos undergo Am^ independent average 
oscillations and R stays more or less constant with Am^ . For the upward going 
neutrinos in addition to average oscillation there is little amount of decay as well 
and hence Y ~ Nup/Ndown is ~ 1- As Am^ decreases to about 0.05 eV'^ the oscil- 
lation wavelength increases - for upward neutrinos it is still average oscillation but 
for the downward neutrinos, the cos(27rL/Ao,5(;) term becomes negative which corre- 
sponds to maximum oscillation effect and the survival probability of these neutrinos 
decreases, and hence R decreases; while the upward neutrinos continue to decay 
and oscillate at the same rate and Y becomes greater than 1. As Am^ decreases 
further, the downward neutrino oscillation wavelength becomes greater than the 
distance traversed and they are converted less and less and thus R increases and 
Y decreases. Below Am^ = 0.001 eV^ the downward neutrinos stop oscillating 
completely while for the upward neutrinos the cos(27rL/Aosc) term goes to 1, and 
R and Y no longer vary with Am^ . 

For the multi-GeV case (fig. 5) the oscillation wavelength is more than the 
sub-GeV case and for Am^ around 0.1 eV^ the cos(27rL/Aosc) term stays close 
to 1 for the downward neutrinos; while the upward neutrinos undergo average 
oscillations and slight decay and Y is less than 1. As Am^ decreases the downward 
neutrinos oscillate even less and the upward neutrinos also start departing from 
average oscillations and hence R increases and Y decreases. Around 0.01 eV'^ 
the downward neutrinos stop oscillation while for upward neutrinos the oscillation 
effect is maximum (A ~ L/2) and the cos{2t^L/\osc) term is ~ -1 and Y stays 
constant with Am^ . As Am^ decreases further the upward neutrino oscillation 
wavelength increases and they oscillate less in number making both R and Y 
approach 1 for Am^ around 0.0001 eV'^. For multi-GeV neutrinos the decay term 
contributes even less as compared to the sub-GeV case. 

We perform a minimization in the three parameters Am^ , sin^ 20 and a. 
The best-fit values that we get are Am^ = O.OOSeF^, sin^ 29 = 1.0 and a = 
0.3 X lO^^eF^. The minimum that we get is 1.11 which is an acceptable fit 
being allowed at 77.46% C.L.. For the 535 days data the best-fit values that we 
get are Am^ = 0.002e\/^ sin^ 2^ = 0.87 and a = 0m23eV^ with a xLn of 4.14 
which is allowed at 24.67% C.L.. Thus compared to the 535 days data, the fit 
improves immensely and the best-fit shifts towards the oscillation limit, the best- 
fit value of the decay constant a being much lower now. It is to be noted however, 
that the best-fit in this model does not come out to be a = 0.0, viz the only 
oscillation limit. In table 2 we give the contributions to from the i?'s and F's 
at the best-fit value of a and for the a = 0.0 case. 

Table 2: The various contributions to the x^min the best- fit value of a and at 
a^O.O 
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Quantity 


a = 0.3 X 10-5 eV^ 


a = 0.0 eV^ 




0.085 


0.021 




0.011 


0.033 




0.48 


0.56 




0.014 


0.073 


e 


0.344 


0.344 




0.176 


0.176 



Thus from the contributions to we see that for the best- fit case there is im- 
provement for the multi-GeV R and Y as compared to the a — 0.0 case. The 
for sub-GeV Y also improves. In fig. 6 we plot the " Xmm with Am^ and 

sin^ 1Q unconstrained. There are two distinct minima in this curve - one for lower 
values and another at higher values of a. The best-fit Am^ in both cases is ~ 0.001 
eV'^ . In this model there are two competing processes - oscillation and decay. For 
lower values of a the decay length is greater than the the oscillation wavelength 
and oscillation dominates. The decay term exp(— aL/i?) is close to 1 and does not 
vary much with the zenith distance L. As a increases the exponential term starts 
varying very sharply with L and the variation is much more sharp for the sub-GeV 
as compared to multi-GeV. This behavior is inconsistent with the data and that is 
why one gets a peak in Ax^ for higher a. As a increases further the ex.-p{—aL/ E) 
term goes to zero for the upward neutrinos and there is complete decay of these 
neutrinos while the downward neutrinos do not decay, the exponential term still 
being 1. Whenever the exponential term is or 1 for the upward neutrinos, the 
wrong energy dependence of this term does not spoil the fit and these scenarios 
can give good fit to the data. Even though fig. 6 shows that the data allows a wide 
range of a, we get the two distinct minima in the Ax^ vs. a curve for high and 
low a values, for both the 535 days (dotted line) and 848 days (solid line) data. 
But while the 848 days data prefers the lower a limit, the 535 days data gives a 
better fit for the high a hmit. The reason behind this is that for the 848 days data 
the R is much higher than for the 535 days data. Hence the 848 days data prefers 
lower a and hence lower suppression. 

In fig. 7 we show the 90% and 99% C.L. allowed parameter region in the Am^ - 
sin^ 29 plane for a range of values of the parameter a. In fig. 8 we show the 90% 
and 99% C.L. contours in the a - sin^ 29 plane fixing Am^ at different values. 
These contours are obtained from the definition < X^min + ^X^: with Ax^ = 
6.25 and 15.5 for the three parameter case for 90% and 99% C.L. respectively. The 
bottom left panel in fig 7 is for the best-fit value of a. For high a (the top left 
panel) no lower limit is obtained on Am^, because even if Am^ becomes so low 
so that there is no oscillation the complete decay of upward neutrinos can explain 



10 



their depletion. As we decrease a the allowed parameter region shrinks and finally 
for a = we get the two parameter limit modulo the small difference in the CL. 
definitions for the two and three parameter cases. The upper right panel of fig. 8 
corresponds to the best-fit value of Am^ . For very low a, even though there is no 
decay, we still have oscillations and that ensures that when Am^ is large enough 
there is no lower bound on a as evident in the fig. 8. For Am^ = lO'^eV"^ the 
neutrinos stop oscillating and hence we get a lower bound on a beyond which the 
depletion in the neutrino flux is not enough to explain the data. 

4 Comparison and Conclusion 

In fig. 9 we show the histogram of the muon event distributions for the sub-GeV 
and multi-GeV data under the assumptions of — oscillation, and the two 
scenarios of neutrino decay for the best-fit values of the parameters both for the 
535 and the 848 days of data. From the figs, it is clearly seen that the scenario 
(a) (big dotted line) (Am^ > O.leV^) does not fit the data well there being too 
much suppression for the sub-GeV upward going neutrinos and too less suppression 
for the multi-GeV upward going neutrinos. The scenario (b) (Am^ unconstrained, 
small dashed line), however, reproduces the event distributions well. However with 
the 848 days data the sub-GeV events are reproduced better as compared to the 
535 days data and the quality of the fit improves. 

The neutrino decay is an interesting idea as it can preferentially suppress the 
upward fiux and can cause some up-down asymmetry in the atmospheric neu- 
trino data. However the intrinsic defect in the decay term exp(— aL/ii^) is that 
one has more decay for lower energy neutrinos than for the higher energy ones. 
Thus neutrino decay by itself fails to reproduce the observed data 0] . If however 
one considers the most general case of neutrinos with non-zero mixing then there 
are three factors which control the situation 

• the decay constant a which determines the decay rate 

• the mixing angle 9 which determines the proportion of neutrinos decaying 
and mixing with the other flavour 

• the Am^ which determines if there are oscillations as well 

If the heavier state decays to a state with which it mixes then Am^ has to be 
> 0.1e\^^ because of bounds coming from K decays [|T^. The best-fit value of a 
that one gets is 0.33 x lO'^eV"^ with the latest SK data. At this value of a the 
^-aL/E |;gj^j2i tends to 1 for the downward going neutrinos signifying that they do 
not decay much. The survival probability goes to (1 — | sin^ 29) which is just the 
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average oscillation probability. In order to suppress this average oscillation the 
best-fit value of sin^ comes out to be small in this picture. For the upward going 
neutrinos, in scenario (a), there will be both decay and average oscillations. If one 
had only average oscillation then the probability would have stayed constant for a 
fixed value of the mixing angle 9. But because of the exponential decay term the 
survival probability drops very sharply as we go towards cos 0=-l.O. The drop and 
hence the decay is more for lower energy neutrinos. As a result the sub-GeV flux 
gets more depleted than the multi-GeV flux, a fact not supported by the data. In 
fact the 848 days data requires the sub-GeV flux to be even less suppressed than 
the multi-GeV flux as compared to the 535 days data and the fit worsens with the 
848 days data. The small mixing signifies that the i/^ has a large fraction of the 
unstable component 1^2 (see eq. (P^). Hence the constant a comes out to be low 
so that the decay rate is less to compensate this. However even at the best-fit a 
of 0.33 xlO~^eV^ the survival probability in the bin with cos0 between -1.0 to 
-0.6 comes out to be 0.15 for E=l GeV, much lower than the value of ~ 0.5 as 
required by the data. Thus scenario (a) fails to explain the upward going neutrino 
data properly because of two main reasons 

• ^ is low in order to suppress the average oscillations of the downward neu- 
trinos 

• the energy dependence of the exponential decay term is in conflict with the 
data 

In the scenario (b), in addition to mixing with z/,-, the unstable component in 
decays to some sterile state with which it does not mix. In this case there is no 

restriction on Am^ and it enters the fit as an independent parameter. We find 

that 

• The best-fit Am^ does not come out naturally to be in the Am^ independent 
average oscillation regime of > 0.1 eV"^, rather it is 0.003e^^. 

• The best-fit value of the decay constant a = 0.3 x 10~^el^^ implying that 
the decay rate is small so that the mixing angle is maximal (sin^ 9 = 0.5). 

• Large values of a giving complete decay of upward neutrinos are also allowed 
with a high C.L. In fact with 535 days data the best-fit was in this region. 

• The best-fit value of the decay constant a is non-zero signifying that a little 
amount of decay combined with Am^ dependent oscillations gives a better 
fit to the data. 
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At the best-fit values of tlie parameters tliere is no oscillation of the downward 
neutrinos so that the cos(27rL/Ao<ic) term goes to 1. The decay term also goes to 1 
signifying that there is not much decay either for the downward neutrinos and the 
survival probability is ~ 1 without requiring the mixing angle to be low. On the 
other hand for the upward neutrinos there are oscillations as well as little amount 
of decay. The sub-GeV upward neutrinos have smaller oscillation wavelength and 
they are close to the average oscillation limit (survival probability ~ 0.5) while 
for the multi-GeV neutrinos the oscillation wavelength is such that one has maxi- 
mum oscillations and the survival probability is less than 0.5. Thus this scenario 
reproduces the correct energy dependence of the suppression - namely sub-GeV 
is suppressed less as compared to multi-GeV neutrinos. The best-fit value of a 
being even smaller now than the scenario (a) the decay term does not vary 

very sharply with the zenith distance L or the energy E so that its wrong energy 
dependence does not spoil the fit. 

The conversion probability of to Vr is given by 

Pu,u^ = ^ sin2 2^{l + exp(-47rL/Arf) - 2 exp(-27rL/Arf) cos(27rL/Ao,c)} (16) 



The value of Pufj,UT integrated over the energy and the zenith angle, for a = 0.3 x 
10~^eV'^ (the best fit for scenario (b)) is 0.33 for sub-GeV and 0.26 for multi- 
GeV. For a = 0.44 x 10~^eV^^ (the second minima in the Ax^ vs. a curve) the 
corresponding numbers are 0.21 and 0.15, while for the only Vp, — Vr oscillation case, 
the corresponding values are 0.37 and 0.26 respectively. The value of Am? and 
sin^ 29 for both the cases is 0.003el^^ and 1.0 respectively. 

The fig. 9 shows that the zenith angle dependence of the scenario (b) is almost 
similar to the case of — Vr oscillation. But the two cases are very different in 
principle. For the oscillation case a larger 9 implies a larger conversion whereas 
in scenario (b) a larger 9 means the fraction of the unstable component is less in 

and the depletion is less. If one compares the conversion probability as given 
by eq. (|l6|) with the one for the z/^ — v-r oscillation case, then the scenario (b) 
considered in this paper would have smaller number of z/^-s in the resultant flux 
at the detector, especially for the larger values of a which are still allowed by the 
data and the two cases might be distinguished when one has enough statistics to 
detect T appearance in Super-Kamiokande |]l9l or from neutral current events |20 



In our paper we have followed the procedure of data fitting as done in []^. 
Thus we use the ratios for which the common systematic errors get canceled out. 
Strictly speaking one should use the absolute number of events and include all the 
correlations between bins and e-like and /i-like events. But the best-fit points and 
the allowed regions are not expected to change significantly We have compared 
the scenarios of neutrino oscillation and decay with the same definition of and 
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for this purpose of comparison neglecting the correlation matrix will not make 
much difference. Apart from the statistical analysis we have given plots of R and 
Y for various values of the parameters. The allowed parameter ranges from these 
plots are consistent with what we get from our statistical analysis. The histograms 
that we have plotted are also independent of our definition of x^- We have checked 
that if we estimate the allowed ranges from the histograms these are consistent 
with what we get from our definition of x^. Thus we agree with the observation 
in ref. that although this method of data fitting is approximate it works well. 
In our analysis we have used only the SK data because it has the highest statistics 
as compared to the earlier atmospheric neutrino experiments. 

Radiative decays of neutrinos are severely constrained 0] and what we consider 
here are the non-radiative decay modes of the neutrino. Models for neutrino decay 



for the scenario (a) are discussed in 0, ^ . For the scenario (b) the unstable state 



decays to a sterile neutrino state and a light scalar. The model described in in 
connection to the solar neutrino problem can be adapted to emulate this scenario. 
A recent paper has discussed how such a model can be constructed. Since 
in scenario (b) the decay products are invisible there are no distinctive signs of 
the decay. Decay of leptons to the light scalar are prohibited from conservation of 
lepton number. Hence it is difficult to constrain these from laboratory experiments 
p3|. Consequences of such a model for astrophysics has been discussed in p4| . 
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Figure Captions 



Fig. 1. The allowed parameter region in Am^ -sin^ 29 plane for the — oscil- 
lation hypothesis for the 848 days data. The solid line is the area allowed at 90% 
C.L. and the dashed line shows the area allowed at 99% C.L. The best-fit point is 
shown. 

Fig. 2. The variation of R and Y with a for the sub-GeV neutrinos (denoted by the 
subscript sg) assuming neutrino decay with Am^ > O.leV^. The curves are drawn 
at fixed values of sin^ ^=0.03 (solid line), sin^ ^^=0.04 (long dashed line), sin^ ^=0.06 
(short dashed line), sin^ ^^=0.08 (dotted line) and sin^^=0.1 (long dashed-dotted 
line). The short dashed-dotted lines give the SK 848 days results within a ±lcr 
band. Also shown are the R and Y at the best-fit point. 
Fig. 3. Same as in fig. 2 but for multi-GeV neutrinos. 

Fig. 4. The variation of R and Y with Am^ for the sub-GeV neutrinos (de- 
noted by the subscript sg) assuming neutrino decay with Am^ unconstrained. In 
these curves the a is fixed at it's best-fit value of 0.3 x lO'^eV"^. The curves are 
drawn at fixed values of sin^ 9—0.7 (dotted line), sin^ ^=0.6 (short dashed fine) and 
sin^^=0.5 (long dashed line). The solid fines give the curves for the best- fit value 
(sin^ 9 = 0.5) of the — Ur oscillation case. The dotted-dashed lines give the SK 
848 days results within a ±lcr band. Also shown are the R and Y at the best-fit 
point. 

Fig. 5. Same as in fig. 4 but for multi-GeV neutrinos. 

Fig. 6. The Ax^ — -)(^ — ■)(^:^^ vs. a with Am^ and sin^ 29 unconstrained for the 
848 days (solid line) and 535 days data (dotted line). 

Fig. 7. The allowed parameter region for the 848 days data in the Am^ -sin^ 29 plane 
for 4 different values of a shown at the top of each panel. The solid and the dashed 
lines correspond to the area allowed at 90% C.L. and 99% C.L. respectively. 
Fig. 8. The allowed parameter region for the 848 days data in the a-sin^ 29 plane 
for 4 different values of Am^ shown at the top of each panel. The solid and the 
dashed lines correspond to the area allowed at 90% C.L. and 99% C.L. respectively. 
Fig. 9. The sub-GeV and multi-GeV ii event distributions vs. zenith angle for the 
various scenarios considered. A^^ is the number of /i events as given by eq. (3) and 
A^o is the corresponding number with survival probability 1. The panels labelled 
SG(535) and MG(535) give the histograms for the sub-GeV and multi-GeV 535 
days data respectively, while the SG(848) and MG(848) give the corresponding 
histograms for the 848 days data. For the both the sets the solid line corresponds 
to the best-fit — oscillation solution, the long dashed line is for the best-fit 
for scenario (a) and the short dashed fine for the best- fit for scenario (b). Also 
shown are the SK event distributions with ilcr error bars for both the sets. 
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